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In relativistic frameworks, given by the Bethe-Salpeter and light-front bound state equations, the 
binding energies of system of three scalar particles interacting by scalar exchange particles are 
calculated. In contrast to two-body systems, the three-body binding energies obtained in these 
two approaches differ significantly from each other: the ladder kernel in light-front dynamics 
underbinds by approximately a factor of two compared to the ladder Bethe-Salpeter equation. 
By taking into account three-body forces in the light-front approach, generated by two exchange 
particles in flight, we find that most of this difference disappears; for small exchange masses, the 
obtained binding energies coincide with each other. 
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1. Introduction 

In this article we sketch out the results of our theoretical study ^ of three-body bound states, 
carried out in full relativistic framework, given by the Bethe-Salpeter (BS) approach [Q] and light- 
front dynamics (LFD) |j3|, 0]. Eventually, the goal is to investigate the bound states of fermions, 
using e.g. a meson-exchange model in the case of 3-nucleon bound states (triton, ^He), and baryons 
as bound states of three (non-perturbatively dressed) quarks interacting via gluons, and calculate 
their electromagnetic form factors. For simplicity however, and as a first step, we restrict ourselves 
to spinless system. For such system, the three-body BS equation with one-boson exchange interac- 
tion was first solved and the results were presented in [|5|]. As far as we know, the LFD three-body 
equation with one-boson exchange kernel has not yet been solved. Previously, the three-body LFD 
equation was solved in ^ 0] for zero-range interaction. Its solution was also found in the relativis- 
tic quantum-mechanical approach for bound and scattering states [|8]], with a phenomenological 
mass operator. 

Though a system of spinless particles is oversimplified, some its properties are rather general 
and their study gives useful physical insight. Especially, the comparison of results found in the BS 
and LFD approaches is fruitful. An important observation is the fact that for two-body system, these 
two approaches give results very close to each other both for the binding energy and also for 
electromagnetic form factors [p^]. The comparison of BS and LFD calculations, for the constituent 
mass m = 1 and for the exchange mass also /i = 1, is shown in fig. |I| The curve NR means non- 
relativistic binding energy (obtained via Schrodinger equation), which considerably differs, even 
for weakly bound system, from almost coinciding with each other BS and LFS binding energies. 
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Figure 1: Binding energy B of system of two spinless particles interacting by one -boson exchange with 
mass /I = 1 v.s. the coupling constant a — / {I6nm^), calculated by BS, LFD and Schrodinger equations. 
The units are set by m = L The figure is taken from [^. 

It is important to note, however, that the actual ladder kernels in the BS approach and in LFD 
are not identical, nor are they given by the same graphs. In the BS ladder kernel, there is no notion 
of time-ordering in the diagrams. On the other hand, the kernel for the LFD equation is given by the 
time-ordered graphs in the light-front (LF) time. Thus, the second iteration in the two-body ladder 
BS equation (the second-order Feynman box graph, figure ^a)), when represented as a set of time- 
ordered graphs, turns into six LF time-ordered graphs, including two so-called "stretched boxes" 
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(a) 




Figure 2: (a) Feynman ladder graph with two exchanges, (b) One of six time-ordered graphs ("stretched 
box"), generated by the ladder Feynman graph (a). 

with two exchange paiticles in the intermediate state. One of such stretched boxes is shown in 
figure ^b). These stretched boxes with two (and more) exchange particles in the intermediate state 
are implicitly included in the BS equation, but they are not generated by iterations of the LF ladder 
kernel; therefore, they are omitted in the corresponding LFD bound state equation. In principle, 
this will cause a difference between the BS and LFD results. However, direct examination of the 
stretched boxes [|ll|, [l2|] shows that they are small. This explains, why the BS and LFD two-body 
binding energies are very close to each other, despite the fact that the kernels are not identical. 

We will see below that in the three-body problem the situation is quite different: even with 
a simple one-boson exchange kernel, the binding energies obtained by BS and LFD significantly 
differ from each other, in drastic contrast to two-body case. This difference will be explained by 
the three-body forces generated by two exchange particles in fight. 

2. Three-body results for ladder kernel 

We consider a scalar field theory, with interaction —g(j)'^(p, and the bound states consisting of 
three particles <p, interacting via the exchange of a 9 boson. Three-body equation is graphically 
shown in fig. ^ 



Figure 3: Three-body bound state equation with ladder kernel. 
The one-boson exchange BS kernel is simply 



K(k,k') 



{k-k'Y-pL^ + ie ■ 



In LFD, this kernel is represented through LF variables Ij^,x (see e.g. eq. {1 .11) from [^). Exactly 
the same kernels were used in the two-body BS and LFD equations. On energy shell, the ex- 
pressions for BS and LFD kernels, being represented in terms of the same variables, are identical. 
Whereas off energy shell they differ from each other. 

The three-body BS amplitude T{p,q;P) (p and q are the Jacobi four-momenta) depends on 
five scalar variables: p^, q^, p ■ q, p ■ P, q P. The square P^ = is fixed by the value of the 
bound mass M. To solve the equation numerically, we transform it to Euclidean space. 
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dimensionless coupling constant a 



Figure 4: Three-body binding energy for one boson exchange with masses pi = 0.01 and pi = 0.5 (units set 
by m = 1). Solid curves are the BS result; dotted curves are the LFD results; dashed curves are the solution 
of the non-relativistic Schrodinger equation. 



To parametrize the three-body LF wave function V''(^12_l,-'':i2;^3_l,-'':3), we introduce transverse 
momenta k± and the ratios x and construct from them LF Jacobi variables: kn^^^n are defined 
for two-body subsystem 12, ^3x,-^3 are momenta of the 3rd particle. The LF wave function also 
depends of five scalar variables: ^12^1 ^3_li ^12± ■^3_Li -^121 ■^^3- 

Thus we have an integral equation, which is four-dimensional in the BS approach and three- 
dimensional in LFD, for the bound state amplitude which is a function of five independent variables 
in both cases. We solve it numerically by straightforward discretization of the internal and external 
variables; we do not make use of any expansion of Y and i//^ in partial waves, nor any other set of 
basis functions. 

For small binding energies = 3m— M^, and for two exchange masses jx = 0.01 and = 0.5, 
our results are shown in figure ^ (the units are set by m = 1). For comparison, we also include the 
non-relativistic result obtained by solving the non-relativistic Faddeev equation in the momentum 
space. Clearly, all three calculations give different results. We emphasize the large difference 
between BS and LFD three-body binding energies, in contrast to coincidence of the two-body 
ones, shown in fig. |l[ 

In this way, new interesting and important features emerge, highlighting differences between 
the BS approach and LFD. In a thiee-body problem the difference between LFD and time-ordered 
iterated BS kernels cannot be reduced to (negligible) stretched boxes: new LFD (non-loop) dia- 
grams appear, which are shown in figure |[ These graphs correspond to three-body forces which 
are not taken into account in the three-body LFD equation with ladder kernel. On the other hand, 
they are implicitly included in the three-body BS equation. 

The three-body potential, resulted from the first-order relativistic corrections due to the first 
eight graphs of figure ^ was first found in Ref . [ 1 3 1 . Its contribution to the triton binding energy was 
calculated in Ref. [14|. However, in Ref. [15] it was shown that this relativistic correction cancels 
with the corresponding relativistic correction to the second iteration of the one-boson exchange. 
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Figure 5: Contribution to three-body forces from two bosons in flight {K\2j, ) and from pair-creation (^i23)- 

Due to this cancellation, sum of these two corrections does not contribute in the Schrodinger equa- 
tion. However, in truly relativistic framework, the full graphs figure ^ (not a first-order relativistic 
correction) should be taken into account. 

The next question is: Can the difference between three-body BS and LFD results be explained 
by the diagrams of figure ^ In order to answer this question, we performed a perturbative calcula- 
tion of the correction to the LFD bound state mass M3 due to these three-body forces. That is, we 
calculated the set of LFD graphs shown in fig. ^ and then found perturbatively the corresponding 
correction to the LFD binding energy. 

The three-body bound state mass squared M| for the exchange masses pL = 0.01, 0.5 and 
1.5 are shown in figure ^. Like in figure Q the solid and dotted curves are our results from the 
three-body BS and LF bound state equations with the two-body one-boson exchange kernel. As in 
figure ^for binding energy, the BS and LFD masses squared differ significantly from each other 
and the difference increases with increase of the binding energy (and, hence, with the coupling 
constant). 

The diamonds with error bars in figure ^ indicate the sum M^^p + AM3, where AM| is the 
correction due to the diagrams of figure ||, calculated perturbatively. The errors indicate the nu- 
merical uncertainty, mainly due to the numerical evaluation of a 12-dimensional integral AM3 ^ 

J \j/^i23 W This correction does indeed shift the LFD three-body mass squared so that it comes 

in reasonably good agreement with the BS results. This is the effect of relativistic three-body forces, 
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Figure 6: Three-body bound state mass squared M| vs. coupling constant a — / (16;rm^) for exchange 
masses pL = 0.01 (top), 0.5 and 1.5 (bottom). The units are set by the constituent mass: m = I. 

dominated here by two exchange-bosons in flight, with a small contribution from the creation (and 
annihilation) of a pair of constituent particle-antiparticles. 

One should of course keep in mind that here we only present a perturbative estimate of the 
correction due to these effective 3-body forces in LFD. Note that the perturbative calculation is 
unexpectedly good even for large correction compared to the value of itself, especially for 
small values jj. of the mass of the exchange particle. The agreement becomes less pronounced as 
the exchange mass increases, once the exchange mass becomes larger than the constituent mass. 
E.g. if the exchange mass is three times the constituent mass, the perturbative correction explains 
about 70% of the difference in M| between the BS approach and LFD. 

Since the corrections are substantial, one could raise the question: How reliable is a perturba- 
tive calculation in this context? To answer this question, one can solve three-body LFD equation 
beyond the perturbative framework, which requires significantly more work. Nevertheless, we have 
ample evidence that the contributions of effective 3-body forces in LFD are substantial in the three- 
body systems, and can explain most of the difference between the results found in the BS approach 
and in LFD. 

3. Conclusion 

We have solved, for the first time, the three-body BS and LFD bound state equations with a 
one-boson exchange kernel and found a significant difference in the corresponding binding ener- 
gies. This difference is absent in the two-body BS and LFD equations [^]. After incorporating, 
in the LFD framework, the three-body forces of relativistic origin, generated by two exchange 
particles in flight, the BS and LFD three-body binding energies become very close to each other. 

In the case of three-nucleon system, there are many other sources of three-body forces. Part of 
them is generated by the transitions N ^ A. In addition, there can be intrinsic three-body forces, 
such as the six-nucleon contact term A^A^A^ NNN, generated by chiral perturbation theory for the 
effective nucleon interaction, or, in QCD, three-body forces due to the triple gluon vertex. In any 
case, the two exchange particles in flight are present and give important contributions. The influ- 
ence of such three-body forces on observables significantly exceeds the experimental precision. 
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